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A general analysis of the consequences of including a fourth neutrino in the standard model matter content, on 
CP violating observables at neutrino oscillation experiments, is presented. Neutrino oscillations in vacuum and 
with matter effects are studied. For the former we update and generalize previous studies on CP asymmetries 
with an additional active neutrino using an updated fit of the PMNS mixing matrix. We study the values of 
the new CP violating phases which maximize the different CP asymmetries in T2K and MINOS-like setups 
aiming to elucidate if the new phases yield measurable effects in the most favorable case. We show that due to 
a combined effect of kinematics and unitarity it is possible to obtain an observable asymmetry in the survival 
channels without violating CPT. For the MINOS-like setup, we find maximum asymmetries in vacuum of the 
order of 2% and 4% for the — > v a and v a — > v T channels respectively. For the T2K-like setup we obtain 
maximum asymmetries of the order of 6% in the survival — > channel. Tree level matter effects enhance 
the former reaching asymmetries of the order of 10% for the — > v e and v c — > v T channels, while the 
fy. — > Vy, survival channel changes slightly depending on the mass hierarchy. Box diagrams with the fourth 
mass eigenstate as a virtual particle were also considered, the corrections to the scattering amplitude being 
negligible. 
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I. INTRODUCTION 

The standard model yields a precise description 
of the fundamental interactions. Nevertheless, in the 
neutrino sector there are still observed phenomena 
whose proper description requires the introduction 
of new elements. Indeed, on one side the tiny neu- 
trino squared mass differences suggests very small 
neutrino masses whose explanation requires the in- 
troduction of new heavy fields; on the other side, 
we have not been able to explain the short-baseline 
anomalies such as the LSND signal fU, the Mini- 
BooNE excesses J2l, and the reactor anomaly JH]. 

The standard picture includes 3 active neutrinos 
whose mixing is described by the Pontecorvo-Maki- 
Nakagawa-Sakata (PMNS) matrix. The mixing an- 
gles and mass differences have been measured in 
different experiments which are designed to maxi- 
mize a desired effect depending on the source, base- 
line and energy. The last parameter to be measured 
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was the mixing angle 6*13 flUt], that was for a long 
time assumed to be zero. 

In this work we explore the consequences on CP 
asymmetries in neutrino oscillations due to the very 
existence of an additional heavy neutrino. As mea- 
sured on LEP, the invisible width of the Z boson 
imposes that an additional active neutrino must be 
heavier than while LEP2 excluded a heavy 

charged lepton, that would be the doublet partner 
of the neutrino considered, up to 100 GeV @]. 
The most attractive scheme would be to consider 
a fourth sequential chiral generation, and indeed, 
there are several motivations to consider this as a 
first attempt to modify the Standard Model (SM). 
These have already been summarized in Q and 
include: new CP Violation (CPV) source for the 
baryon asymmetry of the universe problem lIsT-llDl. 
new perspectives in the Higgs naturalness problem 
or in the fermion mass hierarchy problem IU2T - U6II . 
However, the recent LHC (ATLAS and CMS) re- 
sults put strong limits on the fourth generation Stan- 
dard Model (SM4) since they have excluded at 95% 
a fourth generation down (6') and up quark (if) with 
masses smaller than my < 611 GeV ltl7l - tl9ll and 
m t ' < 557 GeV 1120 - 42311 . Since these quarks cou- 
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pie to the Higgs boson with a strength proportional 
to its mass, they do not decouple from the produc- 
tion of the Higgs boson. The existence of these 
extra fermions, regardless of their mass, would im- 
ply a Higgs boson Mg > 600 GeV. Also precision 
electroweak observables constrain the difference of 
fourth generation quark and lepton masses ll24l - l26ll . 

In this paper though, we will make a discussion of 
the neutrino sector considering a fourth active neu- 
trino regardless of the model in which it is embed- 
ded. Our goal is to test if the two extra CP violating 
phases that appear in the 4x4 PMNS matrix can 
yield a sizable asymmetry in neutrino oscillations, 
hence we scan the whole parameter space for these 
phases obtaining those values that yields the maxi- 
mum CP asymmetry in the different channels. For 
clarity and simplicity, we will start with the results 
obtained for neutrinos propagating in vacuum. We 
point how the somewhat surprising result of hav- 
ing a measurable asymmetry in the survival channel 
v n is obtained due to unitarity and kinemat- 
ical constraints. Then we will proceed in an anal- 
ogous analysis to the case of neutrinos propagating 
through matter for distances and energies that are 
comparable to current experiments. Finally, we will 
calculate the order of the corrections that could be 
induced in the scattering amplitude due to virtual ef- 
fects of the heavy neutrino and its mixing with the 
light ones. 

Our paper is organized as follows: in the next sec- 
tion we discuss neutrino oscillations in vacuum, the 
kinematics of the heavy neutrino and its impact on 
neutrino CP asymmetries. In section III we calcu- 
late the same observables in matter. In Section IV 
we calculate the box diagrams contributing at the 
next order in perturbation theory. Our conclusions 
are given in Section V and we give some details of 
the calculations in an appendix. 



H. NEUTRINOS IN VACUUM 

The neutrino of flavour a, v a , is by definition the 
one that is produced in the weak interactions with 
its charged lepton partner W + — > + v a , where 
a = e, [i o t. These are the interaction eigenstates 
which are related to the propagating states Vi, i = 



1,2,3 through the PMNS matrix U ai as 

K) =XXik>- (!) 

i 

The lagrangian for charged weak currents is given 
in terms of the propagating neutrino state as 

C-cc = ^hal^UaiVLiW^ + fl.C, (2) 

Several considerations are made in order to describe 
a neutrino state at a distance L from the production 
point and at a selected neutrino energy E. First, it 
is assumed that the neutrino mass eigenstate prop- 
agates as a plane wave, then we take the z— axis 
along the neutrino direction and consider that the 
propagating states are ultra-relativistic, i.e., £!j ^> 
rrii. Under these assumptions one can easily find 
the neutrino state after a time t during which it trav- 
els a distance z J^-HH 

K(t)) = U al e-^\ Vl )+U a2 e-^\v 2 ) 

+ ^ 3 e- # > 3 ), (3) 

where = fa. The probability for oscilla- 

tion of initial flavour a to final flavour /3 is given 
by |(i/^|^ Q (£))| 2 and a straightforward calculation 
yields 

P a p = 2TZ{U al U* pi U* a2 Up 2 [e-^ 1 - M - 1]) 
+ 2K(U al U; i U: 3 Up 3 {e-^-^'> - 1]) 
+ 2K(U a2 U; 2 U* a3 Uf] 3 {e- 1 ^-^ - 1])(4) 

where the relation \a + b + c\ 2 = |a| 2 + \b\ 2 + \c\ 2 + 
2K(ab* + ac* + be*) and the unitarity of the PMNS 
matrix are used. Equation <j4j can also be written in 
the compact form 




- 2Y,l{U a iU* pi U* aj U N ) sin (<j> 3 - . (5) 

It is clear that the second term in Eq. (O in general 
will yield a difference in P a p with respect to P a p 
and this effect is driven by the phases in the PMNS 
matrix. 
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A. The PMNS4 x 4 Matrix and the Input Parameters. 

A unitary n g x n g matrix is parametrized by n 2 g 
parameters, out of which 2n g — 1 phases may be re- 
absorbed by rephasing the fields. On the other hand, 
an orthogonal matrix of the same dimension can be 
parametrized by n g (n g — l)/2 angles. So a unitary 
complex n x n matrix will have \ (n g — l)(n g — 2) 
physical phases. For n g = 4 we have 6 rotation 
angles and 3 phases. The PMNS„ x „ matrix can 
be written using the parametrization proposed in 
MM as 

" = Rn g -l,n g Rn g -2,n g ■ ■ ■ Rl,n g ' ... • Rk-l,k 

. . . R23R13R12, (6) 

where is a complex rotation matrix on the ij 
axis and R^ is the rotation without the phase. In 

the case n g = 4 

U = ( ^34(034) x iu 2 4(024, ^3) x w 14 ,(0u, ^2)) 
• (^23(^23) x wi 3 (9 13 ,5) x wi 2 (0 12 )). (7) 

In the calculation of the CP asymmetries below we 
will use for the conventional angles and differences 
of squared masses the best fit points from D"2ll . Al- 
though there are more recent reports on the value 
of 0i3, i.e. ifl, Bl l33i l34l]. th e measured values are 
consistent with those of 0211 and we prefer to work 
with the set of parameters that were used in the fit 
in a single analysis that considers the three neutri- 
nos. As can be seen from Fig. 6 in [32], in this long 
baseline neutrino experiment is not possible to dis- 
entangle the values of 6 and #13, so we decided to 
keep the best fit point which yields (5 = 0. 

The upper bounds on the additional parame- 
ters due to the existence of a fourth neutrino are 
obtained from deviations of the unitarity of the 
PMNS matrix. In this concern, considering elec- 
troweak decays, such as, W decay, invisible Z de- 
cay, test of universality and rare decays, which 
lead to upper bounds for the product (NN^) a p = 
(HV(HVy) a i3, where N is the non-unitary PMNS 
matrix composed of a hermitian (H) and a unitary 
(V) matrix, robust bounds for deviations of uni- 
tarity of the PMNS matrix were obtained in 0~5ll 
which updates previous studies in l36l l37il . We 
use the results reported in 0~5ll to extract the upper 



Parameter 


Value 


Reference 


012 


34.4° 


\M 


013 


9.68° 


[32] 


013 


45° 


[32] 


s* 


0° 


[32] 


Amai 


7.6 x l(T b eV^ 


[32J 


Am| 2 


2.4 x 1(T 3 eV 2 


[32J 


014 


< 3.62° 




0-24 


< 2.29° 




#34 


< 4.21° 




7TI4 


^100 GeV 





TABLE I. Input parameters assuming the Normal Hierar- 
chy (NH). Two variables marked with * are correlated and 
we take their best fit values. 



bounds on the the additional angles due to the exis- 
tence of a fourth neutrino. Relating these bounds 
to the 3 new angles, we obtained: #14 < 3.62°, 
6» 24 < 2.29° and 6> 34 < 4.21°. Notice that we have 
left the extra CP phases, ip 2 and ^3, as free param- 
eters and that from now on we will have a com- 
plex PMNS matrix. We use the input parameters 
shown in Table [I] with Amfj = mf — rrij. For the 
Normal Hierarchy (NH) Am| x = Am 32 + Am^, 
whereas for the inverted hierarchy (IH) we used in- 
stead Ato^ = -2.4 x 1CT 3 eV 2 and 6 13 = 10.99°. 



B. Neutrino Oscillations with an Extra Heavy 
Neutrino 

If a fourth neutrino exists, it is necessarily heavy 
with a mass well above the present energies of neu- 
trino beams, thus the oscillation of a light neutrino 
(y e , u^, v T ) into a heavy neutrino denoted hereafter 
as ve is kinematically forbidden. In spite of this, 
the very existence of such neutrino implies the ap- 
pearance of new phases in the PMNS matrix which 
manifest in observable CP violation effects yield- 
ing asymmetries in the probabilities for neutrino os- 
cillations with respect to those of the antineutrinos. 
Furthermore, as we will show below, the interfer- 
ence with light neutrinos produces interesting ef- 
fects and even in the diagonal channels (surviving 
probabilities) there can be an observable asymme- 
try due to the combined effects of the new phases, 
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unitarity and the kinematics of the oscillations. 

We consider first the propagation of neutrinos. In 
general, the proper calculation of the neutrino state 
at a given time requires to solve the correspond- 
ing wave equation whose complete form is still un- 
known due to the lack of information on the nature 
of the mass term. However, as far as the neutri- 
nos do not interact with other fields while traveling, 
Lorentz covariance allows to write the neutrino state 



as 



ip{t) = ip(0)e- ip ' x = V(0)e 



-i(Et-\p\L) 



(8) 



and the specific form of the state ip(0) is not re- 
quired beyond the fact that it coincides with the 
states produced in weak interactions. The weak 
eigenstates are linear combinations of the neutrino 
mass eigenstates which satisfy the eigenvalue equa- 
tion H\v a ) — E a \v a ), a — 1,2, 3, 4, thus they have 
a component along the fourth neutrino 

\v a ) = U al \l>i) + U a 2\v 2 ) + U a3 \p 3 ) + U ai \Vi) 

(9) 

The fourth component is rapidly damped by the 
kinematics during the propagation. Indeed, in con- 
trast to the light neutrinos, a heavy one propagates 
non-relativistically and we can estimate the size of 
this component at a time t by a simple quantum me- 
chanical calculation. The kinematical suppression 
of this component during the propagations can be 
seen as an imaginary momentum according to 



V to 2 — E 2 . 



(10) 



Performing a non-relativistic expansion it is 
straightforward to show that 



i(Et- \p\L) = -ml 



1 



oS E - 

m 



(11) 



We conclude that for the fourth heavy neutrino in- 
stead of a phase we have an exponential damping 
factor e~ mL which kills this component after a dis- 
tance of order 1 /m. For instance, considering a 100 
GeV mass eigenstate this component disappears af- 
ter a distance of the order of 10 -3 fm. The prop- 
agation effects of a fourth heavy neutrino seem to 
be completely irrelevant. We remark however that 
due to the propagation of and the mixing of the 



neutrinos, at a given time the state \v(t)} contains a 
non-vanishing component of the heavy weak eigen- 
state \ve) and there is a non-vanishing probability 
for the oscillation to this state. Of course, at the 
present beam energies this oscillation is forbidden 
by the kinematics, but then the effects of this non- 
vanishing probability manifest via unitarity. 

In order to explore the possible effects we now 
calculate the oscillation probabilities in the presence 
of a fourth neutrino. The generalization of Eq. (|4} 
is straightforward 



P a p = 2Y,nU al U* 0l U a] Up 3 [e-^-^ - X/P) 
+ Y, 2 nU m U* pz U a4 U 04 (e-^e- mL - 1)). 



Neglecting the term containing the damping factor 
we get the generalization of Eq. (0 as 



P, 



a/3 



i<j 



2 ]T l{U m U* pl U* aj U P] ) sin - ^) (13) 



i<3 



K(U m U^U* A U 0i ). 



There are three major modifications due to the exis- 
tence of a fourth neutrino: i) the new phases modify 
the second term in Eq. (fT~3T > yielding new sources for 
CP violation and CP asymmetries; ii) the appear- 
ance of the last term in this equation; and iii) the 
modification of unitarity relations. In the following 
we explore the consequences of these modifications 
for the asymmetries in neutrino oscillations. 



C. Symmetry Transformations and the Oscillation 
Probabilities 



The symmetry transformations T, CP y CPT 
map the following neutrino oscillations amplitudes 



T 
CP 
CPT 



vp 



Vp -> v a 
D a -> up 
Dp D a , 



thus CPT symmetry requires 

Ppa = Pah 



(14) 
(15) 
(16) 

(17) 
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On the other hand, exchanging v a <-> vp in Eq. (O, 
if the PMNS matrix is complex we get 

P a p + P 0a - (18) 

This means that the weak interactions are not invari- 
ant under the CP (or T) transformation. Combin- 
ing the results in Eqs. ( |17|18t , for a complex PMNS 
matrix we get 

P*f> + P s p, (19) 

and in general it is possible to test CP symmetry 
in weak interactions by considering asymmetries in 
neutrino oscillations. 

Concerning the surviving probabilities, CPT 
symmetry requires 

P(XOt Pqlol (20) 

hence CP asymmetries should not show in the 
survival probabilities. However, the transitions to 
the heavy state are forbidden by kinematics and 
this manifests in observable oscillation asymmetries 
even in these channels. Indeed, probability conser- 
vation requires 

P aa = 1 - Y ( 21 > 

p^a 

where (3 = e, fi, r, E. Similarly 

Paa = 1 - Y P ^ C 22 ) 

P^a 

Thus, using both unitarity and CPT symmetry we 
get 

Y p«e = Y p *e ( 23 > 

P^a p^a 

i.e., CP violation effects cancel in the sum over all 
the channels. However, in spite of being allowed 
by the dynamics, for present beam energies well 
below the heavy neutrino mass, the transition of a 
light neutrino into the heavy one is forbidden by the 
kinematics. On the other hand for a complex PMNS 
matrix from Eq. (foi l we get 



then the combined effect of kinematics, unitarity 
and CP violation is a measurable asymmetry in the 
survival channels even if CPT is a good symmetry, 
e.g. for the muon neutrinos, due to the kinemat- 
ically forbidden transition with probabilities P a E 
and P & e, in an experiment in general we will ob- 
tain 

1 - P„ e - P MT ^ 1 - P m - P m . (25) 

D. CP Asymmetries 

In this section we study the size of the effects of 
a heavy neutrino in neutrino oscillations asymme- 
tries as a function of the involved parameters. The 
sizeable effects of a heavy neutrino in CP asymme- 
tries were shown in ll38ll in the case where the extra 
phases are fixed to 7r/2. This part of our calculation 
updates the input and generalizes results in ll38ll to 
scan the whole parameter space for the new phases 
and to keep the PMNS matrix elements correspond- 
ing to the fourth neutrino independent and within 
the maximal values allowed by unitarity. 

We start by fixing the beam energy and baseline 
to the MINOS (M) and T2K (T) type setup, leaving 
the new phases as free parameters and finding the 
values maximizing the effects. For E we use the 
mean beam energy in both cases. The specific val- 
ues used in the calculation are given in Table [XT] For 
each channel we obtain the values of the new phases 
maximizing the CP asymmetry 

A CP = ^ r , where P± = P af} ± P &p . (26) 

The explicit expressions for the asymmetry in terms 
of the angles and phases characterizing the PMNS 
matrix are lengthy, so we only give in the appendix 
a simplified expression in the channel — > v e , 
where the effects of the new phases are visible. We 
will not consider at all times the last term in Eq. jl3\ 
since the modifications to the asymmetries due to it 
are negligible. 

Our results for the CP asymmetry in the — > v e 
oscillation channel are shown in Fig. [T]and those of 
the — > survival channel are depicted in Fig. 
ED 



PaE 7^ PaE 



(24) 
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Experiment L [km] (E v ) [GeV] 
T2K (T) 295 06 
MINOS (M) 735 3.0 

TABLE II. Values used in the calculation of CP asymme- 
tries 




FIG. 1. CP asymmetry in the v^ — > v e channel for the 
T and M parameters, considering the normal (NH) and 
inverted hierarchies (IH). 



In these plots we can see the maximum values 
reached by the corresponding CP asymmetry and its 
dependence on the new phases. For the — > v e 
channel there are regions in the parameter space 
yielding similar values for the maximum CP asym- 
metries as shown with a line in the upper left plot 
in Figure Q] Also, it is remarkable that the survival 
v ii v ^ channel has the biggest asymmetry when 
evaluated for the T parameters. Comparing between 
the two hierarchies we find a maximum asymmetry 
4 times bigger for the normal hierarchy compared 
with the inverted hierarchy in the — > v e channel 
and for the T parameters. For the remaining chan- 
nels we have a small variation for both, T and M pa- 
rameters. We have studied all the channels finding 
the maximum asymmetries and the corresponding 
phases for each channel given in Table [Till We find 
that the most favoured channels are the — > 




FIG. 2. CP asymmetry in the survival channel t> M — > f M 
for the T and M parameters, considering normal (NH) and 
inverted (IH) hierarchies. 



survival channel for T-like setup and the v e — > v T 
channel for the M-like setup which yield maximum 
asymmetries of the order of 5%. 



T-like 





V» — > Ve 




V e — > V r 


V^i — > 






-93.66° 


-108.77° 


92.81° 


100.30° 


NH 




175.21° 


-70.07° 


-88.26° 


101.39° 




A C p 


0.3% 


0.1% 


0.5% 


5.9% 






-35.09° 


-71.42° 


92.19° 


102.71° 


IH 




-126.11° 


70.35° 


-88.06° 


75.80° 




A C p 


0.1% 


0.1% 


0.7% 


5.7% 
















M-like 










Ufi -> V e 


-¥ V T 


V R — > V'j- 


Vfx — > V[i 






-35.08° 


-110.42° 


92. 30° 


7.85° 


NH 




-126.09° 


-89.35° 


81.77° 


91.52° 




A C p 


2.2% 


0.7% 


4.3% 


0.6% 






79.59° 


-65.38° 


91.80° 


8.68° 


IH 




170.75° 


89.49° 


-79.97° 


-88.44° 




A C p 


2.0% 


0.7% 


3.7% 


0.6% 



TABLE III. Values of the CP asymmetry and phases 992 
and if 3 that maximize it in vacuum. 

For given values of the new phases we get very 
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different results for these setups and this lead us to 
study the asymmetries as a function of the baseline 
and of the neutrino beam energy. We are interested 
in the maximum asymmetry provided by the new 
phases, thus for a given value of E and L we scan 
the phases parameter space keeping those yielding 
the maximum asymmetry. 

Our results for the maximum asymmetry as a 
function of L and E and floating phases cp2, ¥>3 are 
given in Figs. [3] and |4] respectively where we con- 
sider normal hierarchy. 



< 



u 
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1.0 

0.8 

0.6 

0.4 
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0.25 
0.20 
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0.00 



\ 

'■' \ 
'■' \ 

^\ 


0.3 GeV 

0.6 GeV 

1.5 GeV 

3 GeV 

12 GeV 
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400 



600 



800 



0.3 GeV 




0.6 GeV 




1.5 GeV 




3 GeV 




12 GeV ; 























200 400 600 

Baseline [Km] 



800 



FIG. 3. Maximum CP asymmetry in the f M — ¥ i/ e and 
—¥ channels as a function of the baseline consider- 
ing NH. The curves correspond to selected beam energies 
from 0.3 GeV to 12 GeV. 



In order to visualize the way the asymmetry be- 
haves as a function of the baseline once the phases 
are fixed, we consider the T-like beam energy and 
evaluate in the phases corresponding to the maxi- 
mum for the T setup. Our results for the — > 
channel with normal hierarchy are given in Fig. |5] 
were a variation of the order of 6% can be seen 
when varying the baseline. 

From these plots we conclude that high energy 
neutrinos may give a big asymmetry for shorter 
baseline experiments in the appearance channels 
Vu — > v e , — > v T and v e — > v T . CP asymmetries 



0.25 




< 



0.00^ 

" 2 4 6 

Energy [GeV] 

FIG. 4. Maximum CP asymmetry as a function of the 
neutrino's energy for a baseline of 295 Km and 735 Km 
considering the NH. 




< 



100 200 300 400 500 600 700 

Baseline [Km] 

FIG. 5. CP Asymmetry for fixed ipi and ips, using the T 
energy of 0.6 GeV as a function of the baseline with the 
NH. 



can be substantial for low energy neutrinos only for 
specific values of the baselines. We studied these 
observables using the inverted hierarchy obtaining 
similar results. It is interesting that sizeable asym- 
metries in the survival channel — > v^, can be ob- 
tained for low energy neutrinos, for example with a 
baseline between 700 and 800 km and E v around 
0.3 GeV, while there is almost no CP asymmetry 
when considering higher energy neutrinos. 
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III. NEUTRINOS THROUGH MATTER 



fi{t)\v.i) that will be driven by the equation [39- 



As originally pointed out in II39h41|1 . there might 
be non negligible effects on neutrino oscillations 
due to the interactions with matter in which the neu- 
trinos propagate. The interaction between neutrinos 
and matter components: protons, neutrons and elec- 
trons, is feeble but a considerable effect on the os- 
cillation amplitude can be obtained by coherent for- 
ward scattering of neutrinos from many particles in 
the material medium. 

At tree level, the neutral current (Z exchange in 
Fig. |6{a)) yields an identical contribution for all 
neutrino flavours and it only produces a shift in the 
energy eigenvalue Ei. This shift does not affect the 
oscillation probabilities, since it appears as an over- 
all phase factor in the amplitude. The charged cur- 
rent breaks this picture since matter contains only 
electrons and at tree level the only process con- 
tributing is v e e elastic scattering via the u-channel 
W exchange diagram shown in Figure |6j5. At the 
present beam energies well below the W mass this 
diagram yields the following effective interaction 

G F _ 

H eff = -y=e^(l - 75)1/1/7^(1 - 75 )e (27) 
G F _ 

= 77= 67^(1 - 75)e^( 1 ~ 75> 

where a Fierz transformation has been performed 
to obtain the result on the second line. With 
N e w 1.5iVyt/cm 3 , considering a constant density 
p ~ 3 g/cm 3 , the final correction to the energy is 
V2G F N e . 



VeMj,V T 



■ dfj(t) _ m ) 




q,e 



(a) 



q,e 




FIG. 6. Tree level diagrams for neutrino interactions with 
particles in matter. 



Neutrino propagation in matter requires to take 
into account these interactions and the time evo- 
lution of an arbitrary state vector \v(t)) = 



dt 2E 

= H jk f k {t) 



fj(t) + V2G F N e U* e3 U ek h(t) 



(28) 



In 14011 an iterative solution was found. The leading 

rOO 



order for / ■ (i = 1, n) are 



fV(t = 0) = 5 i 



(29) 



Then, assembling the row vectors into anxn matrix 
X that satisfies the equation 



.dX 
'~dt 



XH, 



(30) 



with the boundary condition X(t = 0) = 1 an an- 
alytical solution for ( 1301 is possible for a constant 
electron density N e . This approximation is valid if 
we restrict to L < 750 km. 

We use the results in l40ll and scan the parameter 
space for the new phases to obtain the maximum CP 
asymmetry in a T-like experiment. Our results are 
shown in Figs. |7][8] The asymmetries in general are 
bigger for the channels that involve v e as expected. 
While the results for the other channels show the 
same dependence on the neutrino energy compared 
to the vacuum case, in this channel the low energy 
neutrinos show an enhancement providing a consid- 
erable asymmetry. In the survival channel a varia- 
tion that depends on the hierarchy is found for the 
low energy neutrinos with respect to the vacuum re- 
sult. 



IV. BOX DIAGRAMS WITH VIRTUAL HEAVY 
NEUTRINO 

Although box diagrams are suppressed with re- 
spect to the tree level contributions, enhancements 
of the matter effects can occur when we consider 
a heavy neutrino 14211 . Fig. [9^ shows a box di- 
agram that has already been calculated in ll42ll . 
along with other one-loop diagrams representing the 
whole 0(amJ/M^,) electroweak radiative correc- 
tions to coherent forward neutrino scattering. Here 
we do not consider the whole set of diagrams, since 
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FIG. 7. CP asymmetry in the (a ) z/ M — > u E and the 
v e — > channels with NH (b) — > z/ M channel for 
the normal (NH) and inverted (IH) hierarchy for the T- 
like experiment parameters considering matter effects. 
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FIG. 8. CP asymmetry in the — > u & channel with low 
energy neutrinos considering matter effects. 



we are interested only in diagrams with a virtual 
heavy neutrino, therefore we focus on the box di- 
agrams. We also consider as a possibility the di- 
agrams shown on Figures [9J5 and [5J;, where we 
can have neutrino flavour changing weak interac- 
tions since the virtual propagating neutrino is a mass 



eigenstate, thus we have the possibility of having 
these processes with different neutrino flavours as 
outer legs. 




A B 



FIG. 9. Box diagrams with a virtual heavy neutrino for 
neutrino propagation through matter. 

The evaluations were performed using dimen- 
sional regularization and the Feynman gauge (£ = 
1). We considered the zero external momentum 
approximation, which is a good approximation for 
neutrino energies below 20 GeV. 

The diagrams with the corresponding Goldstone 
bosons turn out to be negligible, since they depend 
on powers of neutrino masses of the first three gen- 
erations. For simplicity we stick to the normal hi- 
erarchy with the mass eigenstate v± being the most 
massive and used the lower bound on m„ 4 = 100 
GeV. We would expect this one to be mostly Ve with 
little components of the three standard flavours, 
since the mixing angles are very small. We used 
for <p2 and ips, that appear in the relevant PMNS 
matrix elements, the values that maximize the CP 
asymmetry of the T parameters for the NH in each 
of the channels. It is only in the case of the flavour 
changing process that the extra CP violating phases 
enter into the expressions. In summary our results 
are the following: 

A) For this diagram the most relevant contribu- 
tion to the amplitude comes from the light neutrinos 
z/j, i = 1,2,3 as intermediate states. For an in- 
coming electronic neutrino the absolute value of the 
amplitude from this diagram is of order 2 x 10 . 
The absolute value for the tree level amplitudes is 
0.4 x 10~ 4 , hence the ratio to the tree level diagram 
is of order 5 x 10~ 3 . Amplitudes with incoming v^, 
v T are even smaller. 

B) This diagram allows for neutrino flavour 
changing processes. For an incoming electronic 
neutrino and flavour conserving diagrams the ratio 
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to the tree level diagram is of the order of 8 x 10 
while for flavour violating diagrams the absolute 
value of the amplitude is of the order of 10 -10 . Con- 
sidering an incoming muonic neutrino with a final 
neutrino either muonic or taonic the order is even 
smaller, namely ~ 10 -12 . 

C) This diagram admits also change in neutrino 
flavour. For an incoming electronic neutrino the ra- 
tio to the tree level process in the ee is ~ 2 x 10~ 3 , 
amplitudes for other flavour conserving processes 
and amplitudes for flavour violating processes are 
at least two orders of magnitude smaller. 

Summarizing this section we conclude that the 
amplitudes of the three box diagrams are at least 3 
orders of magnitude smaller than the tree level am- 
plitude and numerically have no impact in the asym- 
metries in neutrino oscillations. 



V. CONCLUSIONS 

The very existence of a fourth neutrino modifies 
the conventional PMNS mixing matrix introducing 
three new angles and two new phases yielding new 
possible sources for CP violation. Since the mixing 
with light neutrinos is small it is desirable to have 
an estimate of the size of the CP violation effects in 
current neutrino experiment. 

In this work we studied the maximum values that 
CP asymmetries in neutrino oscillations can reach 
when a fourth neutrino is considered, both in vac- 
uum and in the presence of matter. We show that an- 
other consequence of the existence of a fourth neu- 
trino is that it is possible to observe asymmetries 
in the surviving channels even if CPT is conserved, 
due to a combined effect of unitarity, kinematics and 
CP violation. 

In the numerics we set the value of the phase ap- 
pearing in the case of three neutrinos to S — 0, so 
that the imaginary part of the PMNS matrix is due 
entirely to the existence of the fourth neutrino. Us- 
ing the upper bound values for the three extra mix- 
ing angles as obtained form the deviations from uni- 
tarity studies IS, we got 6 U < 3.62°, 6» 24 < 2.29° 
and 834 < 4.21°. We study the maximum asym- 
metries that can be reach in neutrino oscillation as 
functions of the neutrino energy, the distance to 
the detector and the two new phases, scanning the 



whole parameter space for the latter. 

Our analysis shows that, under the above assump- 
tions, in vacuum, the maximum asymmetries that 
can be reached in a T2K-like setup (L = 295 
km, E v = 0.6 GeV) are of the order of of 6% and 
appear in the survival — > channel. Asym- 
metries in other channels are at least one order of 
magnitude smaller. As for the MINOS -like setup 
(L = 735 km, E v — 3 GeV), the maximum asym- 
metries can be reached in the v e — > v T and — > v e 
channel and are of the order of 4% and 2% respec- 
tively. Importantly, for the same baseline but a dif- 
ferent energy within the reach of MINOS, namely 
E v = 1.4 GeV, we get a maximum asymmetry of 
the order of 6% in the surviving — > channel. 

Effects of matter enhance the maximal asymme- 
try in the channels involving electronic neutrinos 
as expected. For the T2K-like setup the maximum 
asymmetry in the — > v e channel raises from 
0.3% in vacuum to 9% in the presence of matter, 
while for v e — > v T it grows from 0.5% in vacuum 
to 10% in a material medium. In these channels re- 
sults for the maximal asymmetries are similar with 
normal and inverted hierarchy. The interaction with 
matter has a small an indirect effect in the survival 
v n v ii channel via unitarity, whose sign depends 
on the hierarchy. It goes from the order of 6% in 
vacuum, independently of the hierarchy, to 8% in 
matter for normal hierarchy, while it decreases to 
4% in the case of an inverted hierarchy. 

The value of the asymmetry in the — > v e chan- 
nel varies with the energy of the neutrino beam in 
general being larger at low energies. For the T2K 
baseline, as shown in Fig. ([8]) the maximum asym- 
metries can be of the order of 70% changing the sign 
when going from beam energies of E u = 0.30 GeV 
to E v = 0.33 GeV. Although not shown explicitly, 
for a fixed energy below 1 GeV we obtain similar 
peaks in the maximum asymmetries for definite val- 
ues of L. For the T2K energy E v = 0.6 GeV the 
peaks appear around L = 550 km. 

We do not report the obtained asymmetries in 
matter in the case of MINOS-like because in this 
case the baseline is at the edge of the validity of the 
approximations used in the solution of the propa- 
gating equations and a deeper analysis beyond the 
scope of this work is necessary in this case. 

The potential contributions to the ve scattering 
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amplitude due to the exchange of virtual heavy neu- 
trinos in box diagrams were found to be negligible. 

Finally we would like to remark that the possibil- 
ity to have an asymmetry in the surviving channel 
Vn — > is a consequence of the very existence of a 
fourth neutrino, thus it would be interesting to try an 
experimental search of this channel since the mea- 
surement of an asymmetry in this channel would be 
a direct proof of the existence of a fourth neutrino 
with important implications in particle physics and 
cosmology. 
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